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Abstract 

The chiral extrapolation of the nucleon mass, M„, is investigated using data coming from 2- 
flavour partially-quenched lattice simulations. A large sample of lattice results from the CP-PACS 
Collaboration is analysed using the leading one-loop corrections, with explicit corrections for finite 
lattice spacing artifacts. The extrapolation is studied using finite-range regularised chiral per- 
turbation theory. The analysis also provides a quantitative estimate of the leading finite volume 
corrections. It is found that the discretisation, finite volume and partial quenching effects can all 
be very well described in this framework, producing an extrapolated value of M„ in agreement with 
experiment. Furthermore, determinations of the low energy constants of the nucleon mass's chiral 
expansion are in agreement with previous methods, but with significantly reduced errors. This 
procedure is also compared with extrapolations based on polynomial forms, where the results are 
less encouraging. 
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I. INTRODUCTION 



There has been great progress in lattice QCD in recent years, associated both with Moore's 
Law and with improved algorithms, which mean that one can work with larger lattice spac- 
ings and still approximate the continuum limit well. The CP-PACS group has devoted 
considerable effort to the study of the masses of the lowest mass baryons and vector mesons. 
This has led, for example, to a comprehensive set of data for the mass of the nucleon in 
partially-quenched QCD (pQQCD), with exceptionally small statistical errors We shall 
exploit this data. 

The remaining barrier to direct comparison with experimental data is the fact that cal- 
culations take much longer as the quark mass approaches the chiral limit. Indeed the time 
for a given calculation scales somewhere in the range to m~^, depending on how hard 

n 

one works to preserve chiral symmetry [2|. As a result there has been considerable interest 
in using chiral perturbation theory (xPT), an effective field theory (EFT) built on the sym- 
metries of QCD, to provide a functional form for hadron properties as a function of quark 
mass . In principle, such a functional form can then be used to extrapolate from the 
large pion masses where lattice data exists to the physical value. Unfortunately, there is 



considerable evidence that the convergence o 



this expansion to be reliable at present |6l-ll3| 



dimensionally regularised xPT is too slow for 



On the other hand, it can be shown that a reformulation of xPT using finite-range regu- 
larisation (FRR) effectively re-sums the chiral expansion, leaving a residual series with much 
better convergence properties [s], Q] . The FRR expansion is mathematically equivalent to 
dimensionally regularised xPT to the finite order one is working s], [l4]. Systematic errors 
associated with the functional form of the regulator are at the fraction of a percent level . 
A tonnal descHpt.on ot the formulafon ot ba^on ,PT u.„g a _un, cutoff (ot FRR) 
have recently been considered by Djukanovic et al. [13]. The price of such an approach is 
a residual dependence on the regulator mass, which governs the manner in which the loop 
integrals vanish as the pion mass grows large. However, if it can be demonstrated that 
reasonable variation of this mass does not significantly change the extrapolated values of 
physical properties, one has made progress. This seems to be the case for the nucleon mass 
11 ] and magnetic moments |l6|, for example, where "reasonable variation" is taken to be 
±20% around the best fit value of the regulator mass. 

In order to test whether the problem is indeed solved in this way one needs a large body 
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of accurate data. This is in fact available for the nucleon, where CP-PACS has carried out 
lattice simulations of pQQCD with a wide range of sea and valence masses. This sector re- 
quires a modified effective field theory, namely partially-quenched chiral perturbation theory 
(pQxPT) [l^ [isl. Formal developments in this field have made significant progress in the 



study of a range of hadronic observables — see Refs. [19|-|2l|, |2J-|27| , for example. 

This large body of pQQCD simulation data is analysed within a framework which incor- 
porates the leading low energy behaviour of partially-quenched EFT. Finite-range regular- 
isation is implemented to evaluate loop integrals, for reasons discussed above. The aim is 
to test whether this approach produces a more satisfactory description of the complete data 
set than the more commonly used, naive extrapolation formulas. 

As we will see, the finite-range regularisation method is able to reproduce the nucleon mass 
with a remarkable level of accuracy using partially-quenched lattice data at only relatively 
large pion mass. Furthermore we are able to determine the low energy constants of the 
chiral expansion of the nucleon mass to a remarkable level of accuracy. For this reason, 
we encourage the generation of partially-quenched data by the lattice community since it 
greatly increases the coverage of parameter space, thus enabling chiral extrapolations to the 
physical point to be performed more accurately. 

This work is a companion paper to 28|] in which we used the same technique to analyse 
the vector meson mass, obtaining the p— meson mass to 1% of its physical value. 

The next section summarises the finite-range regularised forms for the self-energy of the 
nucleon in the case of pQQCD. Section IIIII discusses the data used from the CP-PACS 
Collaboration [l| . We then give details of the chiral fits in Sec. IIV[ Finally, Sec. |V] reports 
the consequences of the fits for the determination of the nucleon mass at the physical point. 



II. SELF-ENERGIES FOR THE PARTIALLY-QUENCHED ANALYSIS 

Theoretical calculations of dynamical-fermion QCD provide an opportunity to explore 
the properties of QCD in an expansive manner. The idea is that the sea quark masses 
(considered in generating the gauge fields of the QCD vacuum) and valence quark masses 
(associated with operators acting on the QCD vacuum) need not match. Such simulation 
results are commonly referred to as partially-quenched calculations. Unlike quenched QCD, 
which connects to full QCD only in the heavy quark limit, pQQCD is not an approximation. 
The chiral coefficients of terms in the chiral expansion (such as the axial couplings of the 
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71 and r]') are the same as in full QCD. Hence, the results of pQQCD provide a theoretical 



extension of QCD [18[. QCD, as realized in nature, is recovered in the limit where the 
valence and sea masses match. 

In this section we explain the form of the finite-range regularised chiral extrapolation 



formula in the case of pQQCD — i.e., the case where the valence and sea c 
necessarily mass degenerate. This work extends on the early work of Ref. 



uarks are not 
1, and mirrors 



our analysis of the vector meson mass in Refs. [28|. 

We restrict our attention to correlators of nucleons containing three degenerate valence 
quarks. However, loop diagrams can (and do) contain baryons which contain non-degenerate 
quarks. For convenience, we introduce the following notation for baryon and pseudoscalar 
meson masses, Mb(/3, Kse^; k^^i, k^j, k^^i) and Mps{(3, K^ea; f^l^^h '^D' where the first two argu- 
ments refer to the gauge coupling and sea quark mass, and the arguments after the semi-colon 
refer to the valence quark values. Throughout the paper it will be convenient to abbreviate 
this by introducing the notation: 

; '^^sea; ^val; ^val; ^valj 

^ = Ais(/3, Ksca! 'tsea; /^val) /^val) 

Mpf = Mp5(/3, Ksea; /tval, ft^val) (1) 
j^non.deg _ ^^^^-J 

) ^sea; '^^sea; ^seaj 

where the superscript unit refers to the unitary data; deg refers to a hadron containing 
degenerate valence quarks; and non-deg refers to the case where the valence quarks are not 
degenerate. 

The derivation of the pQQCD chiral expansion can be described by diagrammatic meth- 



ods 



21| , where the role of sea quark loops in the creation of pseudoscalar meson dressings of 



the nucleon is easily observed. The self-energy considered below (S^v) is the total contribu- 
tion from those pion loops which give rise to the leading non-analytic (LNA) and next-to- 
leading non-analytic (NLNA) terms proportional to F and D in the self-energy of the baryon, 
and also the contributions that arise from the t]' diagrams. Explicitly we write the processes 
as iV -> A^TT -> iV, ^ Att ^ A^, Nr]' N and N Ar]' ^ N. In the limit 

of full QCD these t]' contributions can be neglected because the r]' is heavy and therefore 
decouples from the low energy EFT. The appearance of the unusual term A^ — > Arj' — > A^ is 
a consequence of the fact that intermediate states in the partially-quenched theory are not 
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guaranteed to be physical. In particular, the "A" here denotes a state of spin-3/2 (but not 
isospin-3/2) which is degenerate with the corresponding A state because the hyperfine gluon 
interaction depends only on the spin of the quark pairs. 

For pQQCD in the heavy baryon limit the nucleon self-energy may be expressed as: 

Sat = al,^ + alf^ + a^,^ + a^;^ (2) 

Explicitly we have: 



327r 



-J(M?|, - Mi'')] (3) 



As we will see (figJS]), aj^j^ and aj^^ are typically negative, whereas aj^j^ and aj^^ are typically 
around zero. The parameters a, /3 & 7 are derived from the standard SU{6) couplings |38| . 
Explicitly we take 



a 



2(F + D)2 



2{3F -oy 

7 = -2D 

r = -{3F + Df + 3{D - Ff (4) 
3 

We use the constants F = 0.51 and D = 0.76, which are determined from fitting semi- 
leptonic decays at tree level- e.g., Ref. {32 1. 



The integrals in eq.(l3]) are defined as 

I{Mps,6M) 



PS) 



71 Jq uj{uj + 5M) 
2 k^u\k)dk 



(5) 



where we have used: 



uj{k) = ^Jk^ + Mj,g (6) 

Here Mps can be Mpg or Mp^"~'^^^. We define this along with values for SM explicitly in 
the individual self-energy terms above. 

We study both a standard dipole form factor, which takes the form 



and a Gaussian form factor 



u{k) 



u{k) 



exp 



(7) 




To account for finite volume artefacts, the self-energy equations are discretised so that 
only those momenta allowed on the lattice appear: 



Att I k'dk 
'o 



d^k 



1 /27r 
V \~a 



with 



k 



2n{i,j, k) 



x,y,z 



aK 



(9) 



(10) 



x,y,z 



with i,j,k G Z. The purpose of the finite-range regulator is to regularise the theory as 
kx, ky, kz tend to infinity. Indeed, once any one of kx, ky or kz is greater than ~ lOA the 
contribution to the integral is negligible and thereby ensuring convergence of the summation. 
Hence, we would like the highest momentum in each direction to be just over lOA. For 
practical calculation, we therefore use the following to calculate the maxima and minima for 
i, j, k: 

'lOA a 



2tt 
lOA a 
2% 



N(x,y,z) 



+ 1, 
- 1, 



where [. . .] denotes the integer part. 

We now have the partially-quenched nucleon mass formula which we will use in sec. IIVI 
to fit the CP-PACS data. 
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Volume 


PS I V 


dr- , [fnil 


'-^CT L J 


1.80 


0.1409 


12^ X 24 


8067+^ 


286+*^ 

KJ » ^KJKJ 


288+3 

KJ * KJ KJ 


1 80 


1430 


12^ X 24 




272+2 

KJ . ^ I 2 


280+'' 


1.80 


0.1445 


12^ X 24 


0.694+2 


0.258+1 

— 4 


0.269+2 

— 


1.80 


0.1464 


12-'^ X 24 


0.5471:1 


0.237+1 


0.24812 


1.95 


0.1375 


16^ X 32 


8045+" 


0.196+1 


2044+1*^ 




1 390 


lg3 X 32 


752+2 


1 85+3 




1.95 


0.1400 


16-'^ X 32 


0.6901^ 


174+2 


0.18121^2 


1.95 


0.1410 


16^ X 32 


0.5821^ 


0.16312 


0.169911^ 


2.10 


0.1357 


24^ X 48 


806+2 

KJ • KJ \J KJ Q 


0.1275+^ 


0.1342+1 


2.10 


0.1367 


24^ X 48 


0.755+2 


1203+* 

KJ t i-^KJ\J ^ 


1254+* 


2.10 


0.1374 


24^ X 48 


0.691+1 


0.1157+1 


1203+^ 


2.10 


0.1382 


24^ X 48 


0.5761i^ 


0.10931^ 


0.11291^ 


2.20 


0.1351 


24^ X 48 


0.799+^ 


0.0997+* 


10503+1"^ 

— 15 


2.20 


0.1358 


243 X 48 


0.75311 


0.096611 


0.101313 


2.20 


0.1363 


243 X 48 


0.7051^ 


0.093611 


0.097813 


2.20 


0.1368 


243 X 48 


0.63211 


0.090611 


0.094912 



TABLE I: 



'he lattice parameters of the CP-PACS simulation used in this data analysis, taken 



from Ref. [l|. The superscript unit refers to the unitary data (i.e., where k^^^ = k"^^^ = Ksca.)- Note 
that the errors reported in this table are obtained with our bootstrap ensembles (see text). 

III. THE CP-PACS NUCLEON DATA 

In Ref. [l| , the CP-PACS collaboration published meson and baryon spectrum data from 
dynamical simulations for mean-field improved Wilson fermions with improved gluons at 
four different (3 values. For each value of /3, ensembles were generated for four values of Ksea 
- giving a total of 16 independent ensembles. Table [T] summarises the lattice parameters 
used and Figure [1] is a graphical representation of the unitary pseudo-scalar masses plotted 
against the lattice spacing 0^3 and we note that (Mp™*)^ is a direct measure of the sea quark 



mass. For each of the sixteen ensembles there are five /tvai values considered Thus there 
are a total of 80 {Mp^ , Mpg) data points available for analysis. 
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FIG. 1: The range of sea quark mass (Mp™*)^ and lattice spacing, ar^, covered by the CP-PACS 
data displayed in Table HI (Mp™*)^ is the pseudoscalar meson mass squared at the unitary point; 
i.e., where Kvai = ^sea- The experimental points for the tt,K and "r/^" mesons are also shown for 
reference. 



The analysis in this paper shares many common features with our analysis of the vector 
meson data in Ref. 128 



|. In particular, we consider two methods of setting the scale: the 
string tension (a), and the hadronic scale (rg). We generate 1000 bootstrap clusters for 
all hadronic masses from a Gaussian distribution with a central value equal to the values 
published in [l| (tables XXI, XXII and XXIII 39|) and a FWHM equal to the published 



28 



leads to our 



error. We use totally uncorrelated data throughout, which we argue in Ref. 
statistical errors being overestimates. The values tq = 0.49 fm and y/a = 440 MeV are used. 

Since the action used in [ij] is mean-field, rather than non-perturbatively improved, it will 
have some residual lattice systematic errors of 0{a). We therefore fit the data assuming 
both 0{a) and O(a^) effects, which we investigate in Sec. IIVI 

The physical volume for the /3 = 1.80, 1.95 and 2.10 ensembles is La ^ 2.5 fm, and the 
P = 2.20 ensemble has a slightly smaller physical volume. The associated finite volume effects 
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are incorporated through evaluating the chiral loops by explicitly summing the discrete pion 
momenta allowed on the lattice as described in eq.dH])- 



IV. FITTING ANALYSIS 

A. Summary of analysis techniques 

The philosophy behind our fitting method remains the same as for our investigation of 



the meson spectrum [28|], i.e. we work in physical units when performing our extrapolations. 
We do this so that data from different ensembles can be combined; something which cannot 
be done for the dimensionless data since they correspond to differing lattice spacings. In 
addition, we expect that this approach will benefit from some cancellation of the systematic 
(and statistical) errors. 

The Adelaide approach to chiral fits describes the variation of hadron mass with quark 
mass by a combination of the self-energy term (in this case, Sat) with "constituent quark" 
terms (i.e. polynomials in the valence quark mass). The former accurately describe the 
chiral behaviour and become negligible as the quark mass becomes heavy. Thus we have 

Mn-Ej, = ao + a2{Mplf + a^^Mplf + a,{Mplf . (11) 

In the complete EFT for the partially-quenched theory, there are contributions to this 
polynomial expansion which measure the displacement from the unitarity point (where 
tT^vai = ^sea)- This freedom could be incorporated by extending the terms polynomial 
in the quark mass by 



^2 



{M';if ^ a,{Mplf + a',5l, (12) 



A"vs ' 



and similarly for the third-order term. Here, the notation = (Mpg)"^ — [Mp^^Y is used 
to describe the mass splitting between the sea and valence quarks. In line with our earlier 



work, Ref. 



28 



we make a model assumption by ignoring the terms proportional to 5^. We 
have checked this assumption by confirming numerically that there is no variation of the 
02 coefficient (obtained with the "individual ensemble fit" approach — see next subsection) 



with Mp™*. We found a similar situation in our earlier work, Ref. 28|]. In any case, any 



dependencies below the level of our statistics on are implicitly contained within the chiral 
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self-energies, Eat. Thereby, we use the lattice data to select a preferential regularisation 
scale (and perhaps regulator) which efficiently interpolates between the partially-quenched 
and unitary points. This method has proven very successful in connecting quenched and 



dynamical QCD results over a range of observables, see Refs. [16|, |22|, |23|, |30| for example. 
In effect, the 6 new parameters (at this order) characterising the non-unitarity dependence 
are modelled by a single parameter A. 



Again, following our vector-meson analysis 28[, we contrast this chirally-motivated ap- 
proach with a naive polynomial fitting function, 

= ao + a^iMp^Y + a^{Mplf + a^{M'';if . (13) 

We divide these fits into two categories, "cubic" and "quadratic" depending on whether or 
not the (Mpg)^ is included. 

In fig. [2] we plot the dimensionful nucleon data, M^, and the subtracted nucleon data, 
Mtv — ^n- In the latter, we use a representative value of A = 600 MeV with the dipole form 
factor and the scale set by tq. 

We plot, in fig. [31 the terms which make up S^v (see eq.([2])) for each of the 80 data points 
under consideration, in order to get a feel for their relative size. We also show the continuum, 
physical values for these terms, noting that in this case, the terms involving rj' vanish (as 
required). The curvature in the terms which make up Stv seen in fig. [3] as Mps — )■ matches, 
by design, the LNA and NLNA chiral contributions (see secjlT])- The polynomial chiral fits 
obviously do not reproduce this chiral behaviour. It is this curvature which means that the 
Adelaide approach more accurately reproduces the experimental value of the nucleon mass 
than the naive polynomial approach, predicting a nucleon mass some 60MeV lower than the 
polynomial approach (see later). 

In the next subsection we individually fit Eqs. (!TT]fc [T3|) to the sixteen ensembles (Table [T]). 
Following this we perform a single global fit of the entire data set using Eqs. f lTT] & [T3|) 
modified appropriately by C(a)-style corrections. 



B. Individual ensemble fits 



In this section we treat the sixteen ensembles separately by fitting Eqs. ( fTTl & [T3|) to the 
five degenerate data points {M'^^ , Mp^) from each of the sixteen ensembles. In the Adelaide 
case, we have used our preferred choices of the dipole form factor with the scale taken from 
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FIG. 2: A plot of the degenerate CP-PACS nucleoli data set. Here the scale is set using rg. Both 
the raw data, M^v, and the subtracted data, M^v — S^v, are shown, with the latter shifted vertically 
by 0.2 GeV for clarity. A dipole form factor with A = 600 MeV was chosen in the latter case. 



tq. The other possibilities (i.e., using the Gaussian form factor and taking the scale from the 
string tension) are discussed in Sec. IIV CI We have chosen A = 600 MeV for these individual 
fits which is very close to what turns out to be our preferred value in Sec. IIV CI (Section 
IIV CI discusses in detail the variation of nucleon mass with A.) The fits considered in this 
section are quadratic in the chiral expansion (i.e., we set ag = in Eqs. (ITT] fc [T3l)) since 
cubic fits for the individual fits have 100% error in the and ag coefficients. 

Table [TTl lists the coefficients for both the Adelaide fits and also the naive fits using this 
approach. As expected the leading Adelaide coefficient is always greater than the correspond- 



ing coefficient from the naive fits (ag'^'^' > 



In virtually all cases the 02 coefficient is 



smaller for the Adelaide fits {a'^'^'^^ < ag"™*^). The 04 coefficients are approximately the same 
for both fits but the error in this coefficient is very large, typically 50%. 

We see only that the coefficient is zero within errors in only one ensemble, indicating its 
presence is needed. 
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o ^ „naive „adel „naive „adel „naive „adel 

H "'sea "o 2 2 4 4 

[GeV] [GeV] [Ge\~^\ [GeV"!] [GeV-^] [GeV-^] 



,80 


0. 


1409 


0, 


.9711 


1, 


.03^ 


~4 

-4 


l.llj 


-12 
-14 


1.091i^ 


-0.301- 


-0.29111 


80 


0. 


1430 


0, 


,981^ 


1, 


.04^ 


-3 
-2 


l.llj 


-11 
-12 


1.08111 


-0.29111 


-0.29111 


80 


0. 


,1445 


0, 


.961^ 


1, 


.03^ 


-3 
-3 


1.21j 


-12 
-11 


1.18111 


-0.37111 


-0.37111 


80 


0. 


1464 


0, 


.931^ 


1, 


.01^ 


-3 
-3 


1.27j 


-12 
-10 


1.2311^ 


-0.42111 


-0.41111 



1.95 


0. 


,1375 


1.0011 


1, 


.05^ 


-4 
-3 


1.08111 


1.0711^ 


-0.251^0 


-0.251^° 


1.95 


0. 


1390 


1.0011 


1, 


.06^ 


-3 
-2 


1.041* 


1.021* 


-0.211^ 


-0.201^ 


1.95 


0. 


1400 


0.9911 


1, 


.05^ 


-2 
-2 


l.lllr 


1.081^ 


-0.261^ 


-0.251« 


1.95 


0. 


1410 


1.0111 


1, 


.07^ 


-2 
-2 


1.081^ 


1.051^ 


-0.241« 


-0.231? 



2.10 


0.1357 


1.0411 


1.08: 


h2 
-2 


1.061^ 


1.051^ 


-0.231^ 


-0.231^ 


2.10 


0.1367 


1.0511 


1.10: 


h2 
-2 


1.011^ 


0.991^ 


-0.191^ 


-0.1911 


2.10 


0.1374 


1.0411 


1.10: 


h2 
-2 


1.031^ 


1.011^ 


-0.191^ 


-0.1911 


2.10 


0.1382 


1.0011 


1.06: 


h2 
-2 


1.131« 


1.101« 


-0.2511 


-0.2511 



2.20 


0.1351 


1.041^ 


1.081^ 


1.011 


1.011 


-0.2111^ 


-0.21111 


2.20 


0.1358 


1.1011 


1.1411 


0.87111 


0.86111 


-0.091^« 


-0.081^" 


2.20 


0.1363 


1.0411 


1.0811 


1.03111 


1.01111 


-0.201* 


-0.191* 


2.20 


0.1368 


1.0111 


1.0611 


1.08111 


1.05111 


-0.231* 


-0.221* 



TABLE II: The coefficients obtained from fitting data against Mpg. We list results for both 
the naive and Adelaide fits ( Eqs. (fT3]) & (fTT|) . respectively) for each of the sixteen ensembles listed 
in Table [H A dipole form factor was employed for the Adelaide fits using A = 600 [MeV] and the 
scale was set from tq. 
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FIG. 3: Here we plot the self-energy contributions Eqs.[3]) versus {M™^"'^^^)'^ for the entire degen- 
erate data set (dashed lines are a guide for the eye only). We use the dipole form factor and choose 
an arbitrary value for the Lambda parameter, A = 1 [GeV]. We also include continuum, physical 
data (the solid curves) for the pion processes (the r] case vanishes in the physical limit). 



Figure m is a representative example of one of these fits. It corresponds to the (/3, i^sea) = 
(2.10, 0.1382) ensemble, which is one of closest to the physical point ( fig. [1]). 

With a view to performing a combined, global fit to all sixteen ensembles, we note from 
Table M that there appears to be no discernible trend with the sea quark mass for any of 
the chiral coefficients. However, there does appear to be a lattice spacing effect. Figures ([5] 
& |6]) plot ao,2 against the lattice spacing from tq, ar^- These plots motivate the following 
continuum extrapolation 

„ „cont I individual „ I ^ A\ 

O0,2 — 00,2 +^0,2 '^ro- U^J 

(Note that in sec. IIV CI we investigate both 0{a) and 0{a?) corrections to the chiral coeffi- 
cients.) The results of the fits corresponding to Eq. (IHj) are listed in Table UTTl 

We note that the errors in the Xq"'^*''*'^""' coefficients are around 25%, whereas they are 
more than 50% for X.^'"^*"'^""' . 
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FIG. 4: A plot of Mn versus M^^ for the ensemble (^,Ksea) = (2.10,0.1382). Included are the 
results of the quadratic naive, Eq. (|13p and the quadratic Adelaide, eq. (|lip fits. The scale is set 
from ro, we use a dipole form factor and our preferred value for A (A = 600 [MeV]). 





„cont. 
"0 

[GeV] 


vindividual 
^0 

[GeV/fm] 


xl/d.o.f. 


„cont. 
"2 

[GeV-1] 


j^individual 

[GeV-Vfm] 
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0.97t^ 
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7/14 


Adelaide-fit 


i.uti 


-0.3711'^ 


8/14 


0.961^ 


0.611 


6/14 



TABLE III: The coefficients obtained from the continuum extrapolation of both the naive and 
Adelaide ao,2 values from Table HIl using Eq. (fH|) . 
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FIG. 5: A continuum extrapolation of the oq coefficient obtained from both the Adelaide and naive 
fits Eq. ([MD. 



C. Global fits 



We now analyse the complete set of 80 (Mps, Mm) data points by globally fitting the 



sixteen ensembles of Table [11 As in Ref. 



the idea is that this will produce a highly 



constrained fit and allow us to determine the higher order coefficients in the chiral expansion 
in Eqs.fllllfc [T3|) . the Adelaide scale parameter, A, and the preferred form factor. 

In order to combine data from different ensembles into a single fit, we use the experience 
gained in sec. II V Bl and in Refs. 28|. This tells us that the data's lattice spacing artifacts, 
which are sizable enough for us to discern, lie in the leading coefficie nt, oq. We have checked 
this conclusion by studying combinations of 0{a) and O(a^) terms 40|] in the 02 and higher 
coefficients, but have found that these fits are unstable. Therefore our global fitting functions 
for the Adelaide case is a modified version of Eq. (fTTj) : 



- E 



N 



(15) 
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FIG. 6: A continuum extrapolation of the 02 coefficient obtained from both the Adelaide and naive 
fits Eq. (Hi 



and the global fit function corresponding to the naive case Eq. f|T3|) is 

= {a, + Xr,a-) + a^{Mplf + a,{MplY + a,{Mplf 



(16) 



We consider both "quadratic" and "cubic" chiral fits and also consider both 0{a) and 
0{a^) lattice spacing effects in the coefficient. Hence, for the global fit analysis the 
maximum number of fit parameters in any one fitting function is five. (We discuss the A 
parameter in the Adelaide case later.) Since our data set contains 80 points, the fit will 
hopefully provide highly constrained fit parameters compared to those from the individual 
fitting method (see Table [Tl| • The scale was set using both the string tension and Sommer 
scale, rg. Finally, for the Adelaide method we study both the dipole, Eq. ([7]), and Gaussian, 
Eq. ([8]) form factors. 

The different choices of fitting procedure are summarised in Table IIVI In this table, the 
entries in each column represent a separate possibility, making making a total of 2^ Adelaide 
and 2^ naive fitting types. 

When performing the Adelaide fits, we have to set the A value, see Eqs. (171l8|) . For 
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Fit 


Chiral 




Lattice 


Approach 


Extrapolation 


term in gq 


Spacing from: 












Cubic 


0(a) 


' u 


Adelaide - Gaussian 










Quadratic 




0" 


Naive 









TABLE IV: The different fit types used in the global analysis. 

numerical reasons, we chose a set of trial A values, rather than let it be a free parameter in 
the fitting procedure. In figs. [7]and[8l we show the x^/d-o.f. versus A for each of the fitting 
choices in the Adelaide case, for the dipole and Gaussian form factors, respectively. This 
allows us to study the quality of the fits as a function of A and to fix the best value of A for 
these fits. 

From figs. [7] and [H we see that the behaviour of versus A is not very dependent on 
either how the lattice spacing effects in the coefficient are modelled or on the order of 
chiral expansion. Note also that the best A value, (i.e., the one which minimises x^) does 
not appear to depend on how the lattice spacing effects are modelled. Changing the form 
factor from dipole to Gaussian results in very similar behaviour, except that the A value 
is simply shifted and the dipole gives a slightly lower x^- The biggest effect on the /d.o.f. 
versus A curves is the choice of whether one uses tq or a to set the scale, with rg clearly 
producing the best fits. 

For each of the Adelaide fit choices we have determined the "best" A value, i.e., the one 
which minimises the x^. These are listed in Table IVl We use these values of A to perform 
the sixteen Adelaide fits in Table IIVI The results of these fits, together with the eight naive 
fits, are listed in Tables IVll and IVnl where the scale is set by ro and a, respectively. 

We summarise the results of Tables |Vl] and IVIII below. 

• Fit approach 

We see that the smallest x^/d.o.f. (indicating the best fitting procedure) is given by 
the Adelaide method which uses a dipole form factor. This has the best ^ jd.o.] . in 
every case (independent of how the chiral extrapolation was truncated, how the lattice 
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400 500 600 700 

A [MeV] 



FIG. 7: A plot of /d.o.f against A for the dipole form factor. The dashed horizontal line 
represents increasing from its minimum value by unity for the rg data (i.e. it represents one 
standard deviation). The intercept of this dashed line with the curves (at A =535 and 626 MeV) 
is used to derive upper and lower bounds for the preferred A value for the dipole case. 



artefacts in were modelled and how the spacing was set). Using the Gaussian form 
factor leads to poorer values which are similar to, or worse than the naive approach. 

Chiral extrapolation 

Errors in the cubic chiral coefficient are large compared to their quadratic counterparts. 
However, the cubic fits always produce a non-zero coefficient and they also lead to 
a smaller value than the corresponding quadratic fits. This indicates the need for 
a cubic chiral term. (Note that quadratic chiral fits were preferred in the mesonic 



case |28|.) 



Treatment of the lattice spacing systematics and the fit coefficients 

We see that the fits with 0{a^) rather than 0{a) lead to a lower in the tq case, 

whereas the reverse is true when the string tension is used to set the scale. This 
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FIG. 8: A plot of x^/d-o.f against A for the Gaussian form factor. 



is another indication tliat tliere are 0{a) systematics present in tiie string tension 



data 
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Setting the scale 

The x^/d.o.f. using tq are significantly lower than the a fits, independent of whether 
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we found 



the Adelaide or naive fits were used. In our study of the vector meson 
this to be true for the Adelaide fits, whereas the naive fits had no preference either 
way for ro or a. 



From Figs. [7] and [8] and Tables IVII and IVIIt we see that the best fit choice is the dipole 
form factor using the cubic chiral expansion with O{o?) effects in the Qq coefficient and using 
ro to set the scale. This choice will be used in the next section to determine the central value 
of our nucleon mass prediction and the spread from the other fitting types will be used to 
define the error. 
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Form Factor 


Chiral 


Lattice 


-^bcst 




Expansion 


Spacing 


[MeV] 


Dipole 


Cubic 






)) 




a 


523 


)) 


Quadratic 


ro 


600 


)) 




(J 


550 


Gaussian 


Cubic 




525 


)) 




(J 


475 


)) 


Quadratic 


ro 


550 


?? 


71 


a 


500 



TABLE V: The best A values for each of the Adelaide fits from figs. [7] and [3 Note that best A 
values do not depend on whether 0{a) or 0{o?') corrections in the oq coefficient are used. 

V. PHYSICAL PREDICTIONS 

In this section we extract a value for the continuum nucleon mass, M^v, in the limit of 
physical, degenerate quark masses. We also list the renormalised coefficients (i.e. low energy 
constants in the chiral expansion of the nucleon mass). All our predictions will obviously 
come from the global fit approach of sec. IIVCI since this is a much more highly constrained 
method than the alternative method of sec. IIVBI We obtain our predictions of Mn by 
setting Mpl = Mp°s~'^^^ = M^^^^ = /i^ in Eqs. (HSl [I6] & El) with /i^ being the physical pion 
mass, which we take to be 138 MeV. We also set M^'^ = Ml'"'-'^^^ and M^'^ = M^"""'^''^ in 
Eq. ([3]). In doing this we see that the 77' contributions to the total self-energy (see Eqs. ([2] 
& [3])) disappear in this continuum, physical required. The only remaining term 

involving M^v and M^, is the a]^^ self-energy term in which we set the mass splitting equal 



to the physical mass splitting of the nucleon and A (i.e., 293 MeV |33|). 

We make a physical prediction for each different fitting method (Table IIVP using the 
coefficients (a™"*, 02, 04 & ag) in Tables IVll and IVlTl for the sixteen Adelaide and eight naive 
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Fit 
Approach 


Form 
Factor 


[GcV] [GcVfni-i] [GcVfm-2 


fl2 Obi 

] [GcV-1] [GcV-3] 


ag 
[GeV-5] 








Cubic chiral extrapolation 


ao contains 0{a) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.081=2 -0.231^ 
1.08+', -0.22+2, 
1.021', -0.271^ 


1.2Q+_l -0.51^ 
1.191^° -0.51', 
1.2911° -0.61', 


0.171° 
0.161° 
0.211° 


41 / 75 
47 / 75 
45 / 75 






Cubic chiral extrapolation 


ao contains 0{a?) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.06011^ - -0.621= 
1.05911^ - -0.601= 
0.9991^? - -0.741= 


1.2111° -0.53111 

1.1911° -o.5ii;i 

1.3011° -0.64i;| 


0.171° 
0.161° 
0.221° 


39 / 75 
44 / 75 
44 / 75 






Quadratic chiral extrapolation 


ao contains 0{a) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.1061^ -0.231^3 
1.1011^ -0.221^3 
1.0561^ -0.2811 


1.031', -0.21011? 
1.031', -0.20711? 
1.071', -0.23011? 


- 


43 / 76 
49 / 76 
49 / 76 






Quadratic chiral extrapolation 


ao contains O{o?) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.0881? - -0.631= 
1.0841? - -0.611= 
1.0341? - -0.751= 


imtl -0.20911? 
1.031^ -0.20611? 
1.071^ -0.23011? 




42 / 76 

47 / 76 

48 / 76 



TABLE VI: The results of the global fit analysis where the scale is set from tq. 

fit cases. For the Adelaide fits, we use the relevant preferred value of A taken from Table El 
We list these predictions for M^v in Table IVlllI 

In Figs. [9] and [10] we present a graphical representation of our study of the A dependence 
of Mjv for both the dipole and Gaussian form factors. For the dipole case, we include the 
acceptable range for the A parameter which is represented by two vertical dashed lines. This 
range is defined from our plots of /d.o.f against A (Figs. [7]& [H]) by increasing from its 
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Fit 
Approach 


Form 
Factor 


\J J. z, 

[GcV] [GcVfm-i] [GcVfni-2 


0,2 

] [GcV-i] 


a4 
[GcV-3] 


ag 
[GcV-5] 








Cubic chiral extrapolation 


ao contains 0{a) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.002t;^ -0.171^2 


1.321« 
1.301« 
1.391^ 


-0.71= 
-0.71= 
-0.81= 


0.2811° 
0.27+;° 

' — 10 

0.3311° 


61 / 75 
66 / 75 

62 / 75 






Cubic chiral extrapolation 


ao contains 0{a^) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


0.98611^ - -0.481^ 
0.988tf, - -O.Utl 
0M7t\l - -OMtl 


1.321^ 
1.301^ 
1.391^ 


-0.71= 
-0.71= 
-0.81= 


0.2911° 
0.281^° 
0.331^° 


66 / 75 
69 / 75 
68 / 75 






Quadratic chiral extrapolation 


ao contains 0{a) 






Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.0361^ -0.191= 
1.0361^ -0.171= 
1.0061? -0.221= 


1.081= 
1.081= 
1.111= 


-0.241= 
-0.231= 
-0.261= 


- 


67 / 76 
71 / 76 
69 / 76 






Quadratic chiral extrapolation 


ao contains 








Adelaide 
Adelaide 
Naive 


dipole 
Gaussian 


1.0201^ - -0.491« 
1.0211^ - -0.461« 
0.9871« - -0.581« 


1.081= 
1.071= 
1.111= 


-0.231= 
-0.231= 
-0.251= 




71 / 76 
75 / 76 
75 / 76 



TABLE VII: The results of the global fit analysis where the scale is set from a. 
minimum by unity. 

We summarise the results of this section which are outhned in Table IVIIII and figures [9] 
and [To] below. 

• The statistical errors in the mass estimates are typically less than 1% for the quadratic 
extrapolations and less than 2% for the cubic extrapolations. 

• We see disagreement between all types of fit when different methods are used to set 
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Estimate Form Mn [GeV] [GeV] 

Approach Factor (Scale from ro) (Scale from a) 



Experimental 



0.939 



Adelaide 
Adelaide 
Naive 



Cubic chiral extrapolation ao contains 0{a) 



dipole 

Gaussian 



0.984115 



0.9731it 



1.046115 



0.950111 



0.93811^ 



0.99211^ 



Adelaide 
Adelaide 
Naive 



Cubic chiral extrapolation ao contains 0{o? 



dipole 

Gaussian 



0.96511^ 
0.95611^ 



1.02311^ 



0.93411^ 



0.92311^ 



0.974111 



Adelaide 
Adelaide 
Naive 



Quadratic chiral extrapolation oq contains 0{a) 



dipole 
Gaussian 



1.0061^ 
0.986+1 



1.0761^ 



0.974+« 



0.9591« 



1.02716 



Quadratic chiral extrapolation oq contains 0{a? 



Adelaide 
Adelaide 
Naive 



dipole 
Gaussian 



0.9881® 



0.9691? 



1.0541? 



0.9581^ 



0.9451^ 



1.0081^ 



TABLE VIII: Estimates of M^r obtained from the global fits. (The errors are statistical only.) Our 
experimental estimate comes from a simple average of the proton and neutron masses. 
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FIG. 9: A plot of Afjv as a function of A from the A(ielaide approach using a dipole form factor. 
Recall that the best A value when the scale is set from tq for the dipole form factor is A = 594 
MeV. The two vertical dashed lines define the range of acceptable A values (535 MeV < A < 626 
MeV) obtained by increasing by unity in Fig. [71 



the scale. When the scale is set from tq the mass predictions are always higher than 
when the scale is set from a. 

We also see that the Mn predictions within one particular method (i.e. the Adelaide 
dipole, Adelaide Gaussian or naive method) have a variation in the results of between 
3% and 5%, with the largest variation in the naive mass predictions. This disagreement 
suggests some instability in the fits, possibly because the lattice systematics are more 
complicated than we have assumed. 

The Adelaide method always produces a mass prediction closer to the physical nucleon 
mass. For the cubic fits the Adelaide mass predictions are very accurate compared to 
their naive counterparts. They are typically within two statistical error bars of the 
experimental mass. 
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FIG. 10: A plot of M^r as a function of A from the Adelaide approach using a Gaussian form factor. 
Recall that the best A value when the scale is set from ro for the Gaussian form factor is A = 525 
MeV. 



The variation of in the region of allowed values of A is very small for each different 
fit. Typically of the order of the other uncertainties. 



As with the results of the meson study j28| we conclude by noting all of these points favour 
the Adelaide approach over the naive method. This suggests that the Adelaide method 
should be the preferred method when performing chiral extrapolations and is a significant 
improvement over the naive method. To give the final value for Mn for both the Adelaide 
method and the naive method, we use our preferred choices: the cubic chiral extrapolation 
with O(a^) corrections in oq; to set the scale; and, for the Adelaide method we use the 
dipole form factor. We quote an error that is based on the spread in the mass predictions 
(for the ro case only). We also (for the Adelaide method) include an estimate of the error 
associated with the A parameter which is taken from the vertical dashed lines in Fig. |9l 
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Hence our final mass estimate for the nucleon is 

^Adelaide ^ 965(15)l«l J^MeV (17) 
^Naive ^ i023(15)t^=*MeV (18) 

where the first error is statistical and the second is taken from the fit procedure. The third 
error in the Adelaide case is due to the A parameter. We have not considered any error that 
may be associated with the determination of tq which we take to be 0.49 fm. We note that 
the Adelaide central value would be la (of the combined in quadrature errors) away from 
experiment by simply rescaling tq upwards from 0.49 fm by around 1%. This corresponds 



exactly with what we found in t 



le p— meson mass case 28|]. However, other results have 



hinted at smaller values of ro [35l437l|. To a good approximation, the Mjv estimates for these 
values of tq near 0.49fm can be obtained by a simple scaling of the values in Eqs. f|T71 fc [T8l) . 
We now turn to the renormalised coefficients, Co,2, which are defined from Eqs.(!2] fc [TTl) . 

Mj, = ao + a2{Mplf + a^^Mplf + al^ + al^ + ... 
= ao + a2{M'';if + a^{Mplf 

+ (a^^)W + {al^r\M';if + c,^A{Mplf 

+ (a^^)W + {alJ'\Mplf + c™(MS)^ln(M?|) + . . . 

^ c, + C2{Mplf + c,{MplY 

+CLNA{Mplf + cmlna{M';IY HMpI) + ... (19) 

Again we have used the fact that the r]' terms, o"^^ and o"]^^, disappear in the continuum, 
physical case, and we note that the crJ^Ni'^'Nd. contributions reproduce the LNA and NLNA 
terms respectively once they are chirally expanded. By expanding ajj^ and aj^^ (from 
eq.(l3])) about Mpg = we obtain the values for (a]^Ar)^°'^^ and (aArA)^°'^^ hsted in Tables HXl 
and 1X1 Note that these values correspond to the fit procedures used in Tables I VII and IVIIi 
The A values used are those in Table |Vl From eqf ll9l) we have 

co,2 = ao,2 + (a^^)(°'^) + (a]VA)^°''^- (20) 

Tables HXl and 1X1 list the values of co,2 for each of the fit procedures. Using the same preferred 
fitting method as in the nucleon mass case (cubic chiral extrapolation with 0{a?) corrections 
in ao with ro to set the scale and using the dipole form factor) we obtain 

Co = 0.930(16)i«GeV C2 = 2.61(10)±l?GeV"\ (21) 

27 



where the first error is statistical and the second is from the fit procedure. The cq value is 
very close to the experimental value of M^v as expected. Furthermore, the C2 value compares 
very favourably with who obtained C2 = 2.80(33) (35) GeV~^ (where we note the statistical 
error in 3| was reported at 2a). Note however that this was with full, 2- flavour QCD (rather 
than pQQCD); the reduced error in our case is due to the larger dataset generically available 



in pQQCD. Other work, 3J], calculates a similar quantity they call Ci from low energy tt — N 
fits. This is related to our C2 via 

ci = -C2/4. (22) 



3J] obtains Ci = — 0.91|;j GeV ^, which, using eq. fl2^ . predicts C2 = 3.61^;* GeV ^ Again, 



this estimate is consistent with our result in eq( l21l) . but with significantly larger errors. Note 



however, that the result from [3j] uses experimental (2+1 flavour) data. Thus it is clear that 
there is significant benefit in simulating pQQCD since it generates a larger dataset which 
results in a corresponding reduction in the errors associated with physical predictions and 
low energy constants. 

VI. CONCLUSIONS 

While the computing resources needed to generate gauge configurations with dynamical 
fermions are very significant, the computation of hadron properties with different valence 
quark masses is relatively low cost for any given sea quark mass. Thus, if one can deal with 
the effects of partial quenching in a controlled manner, this approach offers a potentially cost 
effective way to increase the statistical precision of the final, physical results. This philosophy 



has already been successfully applied to the mass of the p meson in earlier work [28[. Here, 
we have estimated the nucleon mass and its corresponding low energy constants from a 
large CP-PACS simulation jl| of partially-quenched baryon masses using the Adelaide chiral 
extrapolation approach. We have shown that this method is a valid approach to the study of 
the chiral properties of the nucleon mass in the partially-quenched theory. As a comparison, 
we have also used a naive (polynomial) fitting procedure, but found that it is a much poorer 
predictor of the experimental number compared with the Adelaide approach. Our predictions 
from both methods are given in Eqs. f|T71 fc [T8|) . 

We have shown that a single global fit of all 80 degenerate nucleon mass points from jlj| at 
different valence and sea quark masses, and at different lattice spacings is possible using the 
Adelaide approach. This fit procedure has only 5 fit parameters and includes chiral effects, 
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Form 
Factor 


[GcV] 


(a]^^)(o) Co 
[GcV] [GeV] 


[GeV-i] 


[GcV-i] 


C2 

[GcV-1] 






Cubic chiral extrapolation 


ao contains 0{a) 






dipole 
Gaussian 


-0.074 
-0.082 


-0.056 0.951:^ 
-0.059 0.94i^ 


1.056 
1.191 


0.339 
0.401 


2.601^ 
2.78t^" 






Cubic chiral extrapolation 


ao contains O(a^) 






dipole 
Gaussian 


-0.074 
-0.082 


-0.056 0.930^;^ 
-0.059 0.918i;^ 


1.056 
1.191 


0.339 
0.401 


2.6111° 
2.7811° 






Quadratic chiral extrapolation 


ao contains 0{a) 






dipole 
Gaussian 


-0.077 
-0.094 


-0.058 0.972t^ 
-0.069 0.9381:^ 


1.066 
1.248 


0.345 
0.433 


2.441^ 
2.711^, 






Vi^ U-dU.! cLLiC LiliicLi t/Atl cLJJUicxtlUil 


Cig CUlllclillD Lyyil J 






dipole 
Gaussian 


-0.077 
-0.094 


-0.058 0.954^8 
-0.069 0.9211? 


1.066 
1.248 


0.345 
0.433 


2.441^ 
2.711^, 



TABLE IX: The renormalised coefficients co,2 using the scale set from tq. 

finite volume and finite lattice spacing effects. The only systematic deviation from nature 
not included is the number of sea quark flavours which is two in the CP-PACS simulation. 
The Adelaide method includes a form factor and an associated scale, A. In this work, we have 
shown that both the type of form factor, and the preferred value of A can be determined. 
As a by-product, we have been able to determine the low energy constant in the chiral 
expansion of the nucleon mass, C2, to a remarkable level of accuracy (see eg. (1211) ). Again 
this emphasises the benefits that can be obtained from the larger dataset that pQQCD data 
affords. 
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Form 






Co 


fo^r ^(2) 




C2 


Factor 


[GcV] 


[GeV] 


[GeV] 


[GeV-1] 


[GeV-i] 


[GcV-1] 



Cubic chiral extrapolation ao contains 0{a) 

dipole -0.051 -0.036 0.9141i= 0.929 0.273 2.521:^ 

Gaussian -0.061 -0.042 0.9001^1 1.078 0.338 2.72i^ 



Cubic chiral extrapolation oq contains C(a^) 

dipole -0.051 -0.036 0.899^;= 0.929 0.273 2.52t| 

Gaussian -0.061 -0.042 0.8861;^ 1.078 0.338 2.72i| 



Quadratic chiral extrapolation oq contains 0{a) 

dipole -0.059 -0.043 0.934t^ 0.977 0.298 2.361^ 

Gaussian -0.071 -0.050 0.9151^ 1.134 0.369 2.58t^ 



Quadratic chiral extrapolation ao contains C(a^) 

dipole -0.059 -0.043 0.918i« 0.977 0.298 2.361^ 

Gaussian -0.071 -0.050 0.900+«„ 1.134 0.369 2.57+? 



TABLE X: The renormalised coefficients co,2 using the scale set from a. 
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